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Abstract 

The goal of this paper is to explore the potential multistabihty of the cli- 
mate of a planet around the habitable zone. We apply our methodology to 
the Earth system, but our investigation has more general relevance. A thor- 
ough investigation of the thermodynamics of the climate system is performed 
for very diverse conditions of energy input and infrared atmosphere opacity. 
Using PlaSim, an Earth-like general circulation model, the solar constant 
S* is modulated between 1160 and 1510 Wm~^ and the CO2 concentration, 
[CO2], from 90 to 2880 ppm. It is observed that in such a parameter range 
the climate is bistable, i.e. there are two coexisting attractors, one char- 
acterised by warm, moist climates (W) and one by completely frozen sea 
surface (Snowball Earth, SB). Linear relationships are found for the two 
transition lines (W^^SB and SB -^W) in {S*,[C02]) between S* and the log- 
arithm of [CO2]. The dynamical and thermodynamical properties - energy 
fluxes, Lorenz energy cycle, Carnot efficiency, material entropy production 
- of the W and SB states are very different: W states are dominated by 
the hydrological cycle and latent heat is prominent in the material entropy 
production; the SB states are eminently dry climates where heat transport 
is realized through sensible heat fluxes and entropy mostly generated by dis- 
sipation of kinetic energy. We also show that the Carnot efficiency regularly 
increases towards each transition between W and SB, with a large decrease 
in each transition. Finally, we propose well-defined empirical functions al- 
lowing for expressing the global non-equilibrium thermodynamical properties 
of the system in terms of either the mean surface temperature or the mean 
planetary emission temperature. This paves the way for the possibility of 
proposing efficient parametrisations of complex non-equilibrium properties 
and of practically deducing fundamental properties of a planetary system 
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from a relatively simple observable. 
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1.1. Planetary atmospheres and extras olar planets 

A very active recent field of research in astrophysical sciences is the ob- 
servational, model-assisted, and theoretical investigation of extra-solar plan- 
etary objects. In less than two decades, the development of new instruments 
have lead to the detection of the first extra-solar planet in the mid 90s to 
more and more refined observations of the characteristics of several hundreds 
of such bodies. The properties of planets vary greatly, in terms of their com- 
positions - gaseous or rocky -, nature of their atmosphere, and of their size 
- ranging in several orders of magnitude. Astronomical and astrophysical 
factors of great relevance include the temperature and the intensity of the 
radiation emitted by the parent star, the orbital parameters, and whether 
or not the planet is tidally locked. A great deal of effort has been directed 
at constraining the combinations of physical configurations potentially com- 
patible with life (habitable zone) , that is to say with the possibility of 
observing water prominently or at least partially in its liquid form at sur- 
face. Obviously, the so-called habitable zone is the setting where we can 
hope to find forms of extraterrestrial life, at least assuming that life nec- 
essarily expresses itself in forms analogous to those we experience on our 
planet. Recently, ESOs HARPS planet finder estimated that just in the 
Milky Way billions of habitable Earth-like rocky planets could orbit around 
faint red dwarfs, with in the order of one hundred in the immediate vicinity 
of our Solar system (Bonfils and coauthors, 2012). We refer the reader to the 



website (|http: / /exoplanet.eu/index.php ) - which is dedicated to collecting 
information on all newly discovered planets and on the related bibliography 



Ferryman, 2011). 



- and recently published books (Dvorak, 2008 Saeger, 2010; Kasting, 2009 



A great deal of interest is now being paid to inferring, studying, and the 
modeling of planetary atmospheres, i.e., in the case of rocky planets, the 
fluid envelope that surrounds the planet and the response to the differential 
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heating due to the absorption of the incoming radiation. Planetary atmo- 
spheres may feature complex chemical compositions and be characterised by 
the relevant presence of phase transitions of some of the components, which 
may have large effects on the radiative properties. The geometry describ- 
ing how the planet and its parent star face each other determines to a first 
approximation the horizontal inhomogeneities, whereas the atmospheric op- 
tical properties, and in particular its opacity to the incoming radiation, has 
an effect on the vertical structure (e.g. Earth's the optically thin atmosphere 
is typically heated from the planets surface). The general circulation of the 
planetary atmosphere is powered, just like a Carnot engine by such inhomo- 
geneities with, on average, heating occurring at higher temperatures being 
then removed in cooler regions. The circulation acts to reduce the temper- 
ature gradients through transport processes, and on average the generation 
of kinetic energy is compensated by its irreversible dissipation through vis- 
cous processes. Irreversible heat transport and dissipation, and processes of 
radiation absorption and emission contribute to the generation of entropy. 
The steady state is realised by balancing energy and entropy fluxes between 
the planetary atmosphere and the surrounding space. Various aspects of this 
closed chain of interdependent processes have been described for terrestrial 



conditions (Lorenz, 1967; Peixoto and Oort 1992 Held and Soden, 2006 



Lucarini 2009 Lucarini et al. 2011), but these arguments could be applied 



to planetary atmospheres more generally. 

The recent discoveries in planetary astronomy and astrophysics are pro- 
gressively affecting classical geophysical sciences (by definition Earth-centered). 
The study of, e.g. Jupiters or Venus general circulation has a rather old his- 
tory, but with the recent and foreseeable future discoveries of exoplanets we 
have the opportunity of facing a vast variety of planetary configurations, so 
that we are going from considering special cases of climates (Earth, Jupiter, 
Saturn, etc), to being able to study a quasi-continuous distribution of cli- 
mates in some parametric space. Following Read (2011), one should note 



that it is possible to reduce the variety of climate settings by adopting the 
fluid-dynamical classical method of similarity, i.e. by defining a set of di- 
mensionless numbers that fully characterise the climate states. When two 
climate states share the same set of dimensionless numbers, they are dynam- 
ically equivalent, so that the statistical properties of one can be mapped into 
those of the other one with simple algebraic operations. At this stage, we 
are left with two additional elements to cope with: 
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68 • how to characterise succinctly a chmatic state, conveying minimal but 

69 comprehensive physical information: the planetary atmospheres are, 

70 in general, turbulent fluids with variability on a very vast of range of 

71 spatial and temporal scales? 

72 • how to verify, at least approximately, the validity of our simulations 

73 and of our theoretical understanding of the planetary circulations we 

74 are able to model? 

75 The two items are closely related to the deflnition of robust observables. 

76 Energy conservation imposes that the incoming radiation onto a planet is 

77 instantaneously equal to the sum of the radiation that is absorbed and scat- 

78 tered by the planet. Assuming steady state, when averaging over a suffi- 

79 ciently long time interval, the radiation absorbed by the planet is equal to 

80 the radiation the planet emits off to space. These quantities allow deflning 

81 the average albedo of the planet and its effective thermodynamic temper- 

82 ature, which constitute the most fundamental description of the properties 

83 of the planet in terms of the flrst law of thermodynamics. The observation 

84 of a planet as represented by a single pixel recorded by a telescope permits 

85 to gather information at this level of detail in terms of its overall macro- 

86 scopic properties. Obviously, the observations provide much more than this, 

87 such as the spectrum emitted by the planet (and by the star). Nonetheless, 

88 in terms of the macroscopic thermodynamic properties of the planet, only 

89 the spectrally integrated quantities are mostly relevant. Under conditions 

90 of steady state, given the inhomogeneity of the incoming radiation and of 

91 the boundary conditions, the vanishing net energy budget at the top of the 

92 atmosphere (TOA) results from the cancellation between regions where the 

93 budget is positive and regions where the budget is negative (low and mid- 
94 high latitudes regions, in the case of the Earth). Steady state argument 

95 imply that the absolute value of the imbalance in either region is equal to 

96 the energy transmitted from the regions where the energy balance at the top 

97 of the atmosphere is positive from those where such balance is negative, and 

98 the transport is performed across the atmospheric fluid envelope through 

99 material transport and horizontal deflection of the incoming radiation (in 
the case of very thick atmospheres). Note that the correct representation of 

101 these large scale properties in state-of-the-art models of the Earth's climate 



100 



102 



is far from being a trivial task (Lucarini and Ragone, 2011) 



103 Moreover, the second law of thermodynamics imposes that the regions 

104 featuring positive TOA energy budget are warmer (i.e. their total emitted 
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radiation is larger) than those where the balance is negative. One can show 
that the entropy production due to such irreversible transport of energy from 
hot to cold regions gives a lower bound to the total material entropy pro- 
duction of the planet. Moreover, it is possible to provide a lower bound to 
the total dissipation of kinetic energy or, equivalently, to the intensity of the 
Lorenz energy cycle (Lucarini et al. , 2011). This discussion implies that if a 



telescope allows for an observational resolution involving more than one pixel, 
able to distinguish between warmer and colder regions, it is possible to de- 
rive fundamental information on the macroscopic thermodynamic properties 
of the planet in terms of the 2nd law of thermodynamics. This brief discus- 
sion suggests that non-equilibrium thermodynamics provides rather powerful 
methods and concepts for analysing the fundamental properties of plane- 
tary atmospheres, for testing models performances, and for deriving bounds 
on the physical properties of the system even when low-resolution data are 



ip,vailable. Following previous work in the field (Pujol 2003; Fraedrich and 



Lunkeit 
2011a|b 



2008 Lucarini 



Li and Chylek 



2009 



2012) 



Lucarini et al. 2010b, 2011[ Pascale et al 



we propose to use this framework in order 
to study planetary atmospheres in a rather general setting of forcings and 
boundary conditions. 

1.2. Habitability conditions and climatic bistability 

As mentioned above, a great deal of interest in the investigation of plan- 
etary atmospheres is directed at studying conditions within or about the 
Habitable zone. Nonetheless, being in the the Habitable zone is a necessary 
but not a sufficient condition for a planet to have hquid water at the surface, 
even if the chemical composition of the atmosphere would in principle allow 
for it. In fact, the paleohistory of our planet provides strong evidence of the 
fact that the astronomical and astrophysical parameters of the Sun-Earth sys- 
tem support two distinct steady states, the warm state (W) characterised by 
widespread liquid water, and one, instead, characterised by the global glacia- 
tion of water and an extremely dry atmosphere, i.e. the so-called Snowball 
state (SB). What we propose here is a thorough investigation of these two 
states in a bidimensional parameteric space, with the goal of detailing the 
region of bistability of the Earth system and analysing the bifurcation points 
and the mechanisms of transitions between the two states. Our investigation 
deals with Earth conditions because this is the only system where an exten- 
sive body of literature observational, theoretical, amd model-assisted - on 
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the Snowball state and the Snowball-Snowfree transitions is available, but 
our scope is much wider. 

Probably the most notable examples of climate change events occurred 
during the Neoproterozoic (period spanning from 1000 million to 540 million 
years ago) , when the Earth is believed to have suffered two of its most severe 



periods of glaciation (Hoffman et al. , 1998) and entered into a SB climate 
state. This period coincided with large carbon dioxide fluctuations, while 
the solar constant (about 1365 Wm-2 in present conditions) is beheved to 
have been 94 % of current levels, rising to 95% by the end of the Neopro- 
terozoic (Gough, 1981; Pierrehumbert et al. , 2011). The two main factors 



effecting concentration of atmospheric CO2 are biotic activity and volcanism. 
Volcanic eruptions bring about very sudden and dramatic increases in CO2 
concentration. This process provides a potential mechanism through which 
the climate state can exit the SB condition, by increasing the opacity of 
the atmosphere and enhancing the greenhouse effect. The biospheric effect 
tends to occur more gradually as the biotic activity and atmospheric com- 
position are coupled so that large fluctuations of the carbon pools take place 
over large time scales. Note that the effect of SB events on the biosphere 
is believed to have been disastrous. Carbon-isotope ratios characteristic of 



Earths mantle (Hoffman et al. , 1998; Kennedy et al. , 1998) rather than of life 



processes recorded immediately below and above the glacial deposits imply 
that oceanic photosynthesis was effectively non-existent during SB events. 
The result of this and anoxic conditions beneath the ice should have lead to 
the disappearance of most kinds of forms of life except bacteria. The final 
disappearance of SB conditions since the Neoproterozoic may have been the 
main contributing factor in the development of complex multi-cellular life 
that began around 565 million years ago. 



Based on the evidence supported by (Hoffman et al. 1998 Hoffman and 



Schrag, 2002), it is therefore expected that the Earth is potentially capable 



of supporting multiple steady states for the same values of some parameters 
such as the solar constant and the concentration of carbon dioxide, which 
directly affect the radiative forcing. It is important therefore to explore this 
hypothesis, due to the relevance for the history of our planet but also to 
help understand other planets capability for supporting life. Initial research 



using simple 0-D models (Budyko, 1969; Sellers, 1969), 1-D models (Ghil 



1976) as well as more recent analyses performed using complex 3-D gen- 



eral circulation models (Marotzke and Botztet 2007 Voigt and Marotzke 



2011 Pierrehumbert et al. , 2011), provide support for the existence of such 
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bistability. The SB— J-W and W— J-SB transitions tend to occur in an abrupt 
rather than a smooth transition. The main mechanism triggering such abrupt 



transitions is the positive ice-albedo feedback (Budyko, 1969; Sellers, 1969) 



Such a feedback is associated with the fact that as temperatures increase, the 
extent of snow and ice cover decreases thus reducing the albedo and there- 
fore increasing the amount of sunlight absorbed by the Earth system. Con- 
versely, a negative fluctuation in the temperature leads to an increase in the 
albedo therefore reinforcing the cooling. The presence of such catastrophic 



climate shifts (Arnol'd, 1992) suggest the existence of a global bifurcation 



in the climate system for certain combinations of its descriptive parameters 



(Fraedrich, 1979). The loss of stability realized in the W— J-SB and SB— J-W 



transitions is related to the catastrophic disappearance of one of the two 
attractors describing the two possible climatic states, as a result of a set of 
complicated bifurcations. 

Starting from present conditions, the most obvious physical parameters 
to modulate in order to bring about the transition to the SB state is the solar 



constant. Even if other model experiments (Voigt and Marotzke, 2011) show 



that the decrease in CO2 alone can bring about transition to the SB state, 
this requires approximately an 80% decrease in CO2 concentration, compared 
to a decrease of less than 10% for the solar luminosity. The Neoproterozoic 
however highlights the importance of considering the changes in the CO2 
levels as a mechanism for the transitions to and from the SB climate state, 
and therefore the dramatic impact it can have on the overall state of the 
climate system. It is therefore interesting to alter both the solar luminosity 
and the atmospheric opacity as these are two important parameters affecting 
the overall properties of the system. If one wants to explore extensively the 
parametric space of climate steady states, it is therefore necessary to consider 
a wide range of values for both of these parameters. 

Using PlaSim, a general circulation model of intermediate complexity 
(Fraedrich et al. , 2005), we study the climate states realised when the solar 
constant is modulated between 1160 Wm~^ and 1510 Wm~^ and the values 
of [CO2] are varied beetwen 90 to 2880 ppm. Our aim here is an empiri- 
cal (i.e. based on model simulations) reconstruction of the global structural 
properties of the climatic attractors. For both W and SB states we com- 
pute surface temperature, material entropy production, meridional energy 



transport, Carnot efficiency (Johnson, 1997) and dissipation of kinetic en- 



ergy and propose empirical relationships in the parametric plane (5'*,[C02]). 
We will look for an empirical relation for the two transition lines (W— J-SB 
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217 and SB— j-W) in the parametric plane between S* and the natural logarithm 

218 of [CO2] which marks the boundaries of the hysteresis in the climate system. 

219 The aforementioned quantities are used to explain changes in large-scale cli- 

220 mate behaviour and the effect of climate change on features such as strati- 

221 fication and baroclinicity in order to understand changes in the meridional 

222 heat transport across the parameter range. It will also be shown that the 

223 Carnot efficiency has a key role in defining the stability of the system which 

224 is related to abrupt climatic shifts. 

225 The paper is structured in the following way: in Section |2] we describe 

226 details of the non-equilibrium thermodynamics of the climate and the diag- 

227 nostic tools used. Section |3] is dedicated to the description of the PlaSim 

228 climate model and of the experimental setup. In Section |4] we discuss the 

229 results of the simulations, in Section |5] we propose parametrisations of the 

230 main non-equilibrium properties of the system as a function of the mean sur- 

231 face temperature or emission temperature, and in Section |6] we present our 

232 conclusions. 



233 2. Non-equilibrium Thermodynamics of the climate 
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In this section we recapitulate some thermodynamic properties of the 
climate system and introduce the notation used throughout this chapter. 



We follow what presented in Lucarini (2009). If the climate system is en- 
compassed by a domain Q, the total energy budget is given by E{Q) = 
P(i7) + K{Q), where K represents the total kinetic energy and P is the 
moist static potential energy, which includes contributions from the thermal 
(including latent heat) and potential energy. The time derivative of K and 
P can be found to be = —D + W and P = + D — W, where D is the 
dissipation and therefore always positive, W is the instantaneous work done 
by the system and \1/ which is the heating due to convergence of turbulent 
heat fluxes and radiative heat, such that E = The dependence on Q 
has been dropped for convenience. As soon as the climate is considered as a 



non-equilibrium steady state system (NESS, see Gallavotti (2006)), we have 

that over long time scales E = P = K = (the bar indicates averaging 
over long time peri ods). Let us d efine Q as the local heating rate so that 



Q = p(e^ — V ■ if) (Lucarini 



2009 ) where e > is the local rate of heating 



due to viscous dissipation o: 
turbulent heat fluxes plus radiative energy fluxes 
into the domain f2+ , where Q = > 0, and Q 



kinetic energy and H is given by the sum of 

We split the domain Q 
in which Q = Q+ < 0. 
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253 Therefore, we find that and Q integrated over f2 equal the derivative of 

254 the total heating due to dissipation, D and the convergence of heat fluxes 



^! + D = P + W 



dVpQ^ + / dVpQ- = $+ + $- = $ (1) 



255 Where the quantities and $^ are positive and negative at all times, 

256 respectively. Since dissipation is positive definite, —K + W = D = P + W = 

257 W = $+ + $^ > 0. On spatial scales far smaller than Q itself, it is practical 



258 to assume local equ i 



DeGroot and Mazur 



ibriu m (local thermodynamic equilibrium hypothesis, 
(1984)) so that locally Q = sT with s the time derivative 



260 of the entropy density. The total rate of change of the entropy of the system 



261 is: 



262 
263 
264 
265 
266 
267 
268 
269 
270 
271 
272 



274 
275 

276 
277 



s 



dVpQ+/T+ / dVpQ-/T 



dVp\s- 



dVp\s- 



(2) 

where S+ > and S < 0. Using equation (2) and assuming that the 
Earth system is in a steady state, over a long time average, S+ = — E~ as 
S = 0. Therefore, 211+ = J^dVp\s\, so that 11+ measures the absolute value 
of the entropy fluctuations throughout the domain. When integrating over 
the whole domain and considering long time average, we have the follow- 
ing equivalent expressions for the thermodynamic quantites: $+ = S+0+ 
and = S^0^ , where G'^ and 0^ are the time and space averaged tem- 
peratures of the r2+ and domains respectively. Since |S+| = |S~| and 
|$"*"| > |$~|, it can be shown that 0+ > 0~, i.e absorption typically occurs 



at higher temperature than release of heat (Peixoto et al. , 1991 Johnson 



1997). The work done by the Car not engine of the climate system is found 
to be, W = ■?7<l>+, where 



V 



(3) 



can be defined as the Carnot efficiency of the system. As shown in Lorenz 



Lorenz (1967) and clarified in Johnson (1997) (2000) the long term average 



of the work performed by the system is equal to the long-term average of the 
generation of available potential energy, as typical of forced-dissipative steady 
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states. The Earth exists in a steady state maintained far from equihbrium 
by net radiative heating in the warm region (low latitude in our planet) and 
net cooling at the cold regions (high latitude in our planet), which are com- 
pensated by large scale transports performed by the planetary atmpsphere. 
This gives rise to ongoing irreversible processes, including phase transitions 
(H2O in the case of our planet) and frictional dissipation, which are charac- 
terized by a positive entropy production. The entropy production due to the 
irreversibility of the processes occurring within the climatic fluid is called the 
material entropy production, Smat and can be written in general terms as: 



Smat= [ '^dV+ I FsH-vl-dV+ I pLH-vl-dV (4) 

where the first, second, and third terms on the RHS are related to the dis- 
sipation of kinetic energy, and to the transport of sensible and latent heat 
respectively. We now wish to link the terms of the entropy budget in eq. ^ 
with those of the entropy production in eq. (|4]). The second law of ther- 
modynamics states that the entropy variation of a system at temperature T 



[Receiving an amount of heat 5Q is larger than or of equal to 5Q/T Landau 



and Lifshitz (1980). In this case: 



In dVpQ) ] ^ + 
J^dVpTj [ (0) 



(0) 



(e++e-) (e++e-) 

2 2 



(5) 



where Smati^) is the long-term average of the material entropy production, 
Smini^) is its lower bound, i.e. the minimal value of the entropy produc- 
tion compatible with the presence of a Lorenz energy cycle with average 
intensity and (0) is the density averaged temperature of the system. 
The approximation holds as long as we can neglect the impact of the cross- 
correlation between the total net heat balance and the average temperature 
and we can assume that (0) can be approximated by the mean of the two 
Carnot temperatures 0"*" and 0~. We can explicitly write Smin{^) as: 

W r/$+ 



'S'mm(^) 



(0+ + 0-)/2 (0+ + 0-)/2 
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©+ ^ ri ^ ^ 

'(e+ + 0-)/2 l-r//2 ' ^ ^ 

304 where the last approximation holds as long as r/ -C 1, which is typically 

305 the case. Therefore, r] sets also the proportionality factor relating the lower 

306 bound to the entropy production of the system Smin{^) due to macroscopi- 

307 cally irreversible processes to the absolute value of the entropy fluctuations 

308 inside the system due to macroscopically reversible heating or cooling pro- 

309 cesses. Note that if the system is isothermal and at equilibrium the internal 

310 entropy production is zero, since ?7 — )• 0. The lower bound to the material 

311 entropy production corresponds to the contribution coming from the dissi- 

312 pation of kinetic energy through viscous processes. Therefore, the average 

313 material entropy production can be expressed as Smat = Smin + Sexc, where 

314 Sexc is the excess of entropy production with respect to the minimum, which 



315 results from the heat transport down the temperature gradient ((Lucarini 



2009)). We can define: 



^^Sex^. JndViFsH_+FLH)-^i^) 

Sm^n iW/{e)) - ^ 

as a parameter of the irreversibility of the system, which is zero if all the 
production of entropy is due to the unavoidable viscous dissipation of the 

mechanical energy. As Smat ~ 77S+(1 + a), we have that the entropy pro- 
duction is maximized if we have a joint optimization of heat transport and 
the production of mechanical work. Note that, if heat transport down the 
temperature gradient is very strong, the efficiency r] is small because the dif- 
ference between the temperatures of the warm and cold reservoirs is greatly 
reduced (the system is almost isothermal), whereas, if the transport is very 
weak, the factor a is small. The parameter a introduced above is related to 



the Bejan number Be a.s Be = a + 1 (Paoletti et al. , 1989). 



327 3. Experimental setup 

328 We study dynamical and thermodynamical properties of the Earths cli- 

329 mate using the PlaSim climate model. This is motivated by the fact that 

330 our study requires altering very extensively some parameters of the climate 

331 system and producing many simulations. Therefore we need a model which 

332 is flexible and fast to run rather then a state-of-the-art model encompassing 
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jfis ma ny processes occurring in the Earth as possible. PlaSim (Fraedrich 



et al. , 2005 ) is a climate model of intermediate complexity, freely avail- 
able at htpp: //www. mi. uni-hamburg.de/plasim. Its dynamical core is for- 
mulated using the primitive equations for vorticity, divergence, temperature 
and the logarithm of surface pressure, solved using the spectral transform 



low, moist (Kuo 1965 



1978 Stephens et al. 



1974) and dry convection, cloud formation (Stephens 



1982 



method (Eliasen et al. , 1970 Orszag, 1970). Unresolved processes for long 
(Sasamori, 1968[ ) and short (Lacis and Hansen, 1974) wave radiation, shal- 



Slingo and Slingo 1991) and large scale precipi- 



tation, latent and sensible heat boundary layer fluxes, horizontal and vertical 



diffusion (Louis, 1979; Louis et al. , 1981 Laursen and Eliasen, 1989) are pa- 



rameterized. The model is coupled to a 50-m deep mixed layer ocean which 
contains a thermodynamic sea-ice model. The advantage of using a slab 
ocean as opposed to a full ocean is that it allows for the climate system to 
reach a stead state in less than 35 years after a change in e.g. the solar con- 
stant. With full ocean coupling, the integration time of the model and the 
time needed to reach a steady state would be an order of magnitude larger 



(Voigt and Marotzke, 2011). We wish to emphasize that whereas most state- 



of-the-art general circulation models feature considerable energy imbalances. 



as highlighted by Lucarini et al. (2011), the energy bias is of the order of 
0.5 Wm~^ in very extreme climate conductions and less than that for more 
usual choices of climatic parameters. Moreover a tested entropy diagnostic 



is available ( Fraedrich and Lunkeit , 2008 ) , thus making it well suited for this 



work. The model is run at T21 resolution (approximately 5.6° x 5.6°) with 
10 vertical levels. Changing S* over a wide range of [CO2] of 90, 180, 270, 
360, 540, 720, 1080, 1440, 2160 and 2880 ppm, we are able to reconstruct the 
SB and W climate states. The procedure occurs as follows for each of the 
considered values of [CO2]: 

1. the model is run to a W steady state for 100 years with S* equal to 
1510 Wm-2; 

2. S* is decreased by a small amount for each value of CO2 and the model 
run is continued until a steady state is reached; 

3. step 2 is repeated until S* is reduced to 1165 Wm~^; the point of 
W— J-SB transition is noted down; 



367 4. the reverse operation is then performed with S* increased step by step. 
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368 up to the value of 1510 Wm ^, each time allowing the system to reach 

369 a steady state; the point of SB— j-W transition is noted down. 

370 Further to this, we identify the position of the transition to a higher resolution 

371 than the rest of the parameter range in the direction of S* . For values of 

372 S* within 10 Wm~^ before the transition, S* is decreased in intervals of 1 

373 Wm~^, each time permitting 50 years for the system to reach a steady state, 

374 until after the transition is observed. 



375 4. Hysteresis, bistability and regime boundaries in a parametric 

376 space 

377 ^.i. Temperature and entropy production 

378 Initially, the focus is put on analysing the parametric plane ([CO2], 5"*), 

379 which in the following shall be referred to as the CS space, as a function 

380 of global mean surface temperature, Tg. The transition zones between the 

381 main climate states are clearly defined from the temperature profile. Note, 

382 the qualitative properties of the climate system in the CS space, namely the 

383 presence of bistability or of just one of the SB/W states, can be reconstructed 

384 from any observable of the climate state, but it is most instructive to select 



385 first the surface temperature because it is also practically the most relevant. 



387 



400 
401 
402 



386 The temperature profile through the CS space is illustrated in Figures 2 (a 



and 2(b)[ We can identify two main climatic regimes, observed as two dis- 



388 tinct manifolds ([CO2], 5**, T^) and characterized by a sharp change in the 

389 profile of Tg when jumping from one manifold to another. We refer to these 

390 as the upper and lower manifolds, representative of the W and SB regimes 

391 respectively. As would be expected, there is a monotonic increase of tem- 



392 perature with increasing CO2 or S* on both manifolds Voigt and Marotzke 



393 (2011); Pierrehumbert (2005). The temperature range on the SB and W 

394 manifolds are 212 K-242 K, and 254 K-326 K respectively, over the paramet- 

395 ric space. We see that the temperature range 242 K-254 K is not permitted 

396 by the climate system. Note, due to the different temperature ranges the 



397 colour scaling of figures |l(a) and |l(b) is a factor of 4 different, with both 



398 scales starting from the same lowest value. 

399 The temperature range of the bistable region in the SB and W regimes are 
218 K-242 K and 254 K-300 K respectively, meaning that the rate of change 
of surface temperature over the same range of S* and [CO2] is approximately 
twice as large in the W regime with respect to the SB regime, and that the 
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surface temperature difference between the two manifolds ranges between 40 
K and 60 K for identical values of S* and [CO2]. The W states (upper man- 
ifold) exists only in the region of the CS space above the W— )-SB transition 
line (the position of this line is expressed as S* = S*^^^) whereas the SB 
states (lower manifold) only in the CS region below the SB— )-W transition 
line (the position of this line is expresses as S* = S^^^jj). Such lines, which 
are well separated and approximately parallel, are illustrated as solid and 
dashed purple dashed lines on figures 1(a) 8(b) and have been found within 



an accuracy of 2 Wm~^ of the solar constant. The solid purple lines indicate 
the active transition, dependant on which manifold the climate system lies in 
at that moment i.e the active transition for the SB and W states are SB— )-W 
and W— >SB respectively. The dashed lines illustrate the location of the inac- 
tive transition, when the climate system exists in the alternative state. The 
bistable region is therefore located between the dashed and solid purple lines. 
As a result, a property of the system is that regardless of which combination 
of [CO2 ] and S* is used, the transition from one state to another always 
occurs at almost exactly the same temperature. This indicates that the cli- 
mate system is "blind" to the mechanism of forcing. The position of the two 
boundaries can be parameterised in terms of S* and CO2 concentration as: 



'-'sbw 



^sbw 

iogio[co2] + a 



sbw ) 



wsb 



,blogio[C02] + C, 



wsb 



(8) 



where asbw ~ ciwsb ~ —70 Wm~^, Csbw ~ 1630 Wm~^ and C^sb ~ 1440 Wm~^ 
for the transition SB— J-W and W— )-SB respectively and [CO2] is expressed in 
ppm. The size of the bistable region, which we define as B, along S* can 
therefore be defined by the difference between Cgbw and C^sb '■ 



B = C. 



sbw 



wsb- 



(9) 



It is found that B is approximately 200 Wm ^. The displacement between 
the position of the boundaries gives a precise measure of the hysteretic prop- 



erties of the climate (Budyko 1969 Sellers 1969 Voigt and Marotzke 2011 



Lucarini et al. 2010b) since it indicates the size of the overlap between the 



two manifolds in the CS plane. The presence of a bistable region implies that 
when we change the values of S* and [CO2] from an initial to a final value, 
the final steady state depends on the initial steady state and on the path of 
change of S* and [CO2]. Let us assume that we start from an initial point 
{Sq, [C02]o) in the bistable region and in the W state. Let us also assume 
we perform a close path of variation of S* and [CO2] so that Sq = 3% and 
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[C02]o = [CO2]/ . If the path does not cross S*^^^ the final state will be 
identical to the initial one, that is, in a time-averaged sense: 



Ts{S;,[C02]o)=TsiS},[C02 



J/J 



(10) 



On the other hand, if the closed path crosses the transition line to the second 
manifold, the final state will be different from the initial: 



T,{S*„[C02]o)^T,iS},[C02 



J/J 



(11) 



If, furthermore, the closed path crosses first 5^^^ and then S**^^, then again 
Ts{Sq, [C02]o) = Ts{S*j-, [CO2]/), since the system has performed first a 
W— J-SB, and then a SB— J-W transition. The same applies starting from a 
SB state and exchanging SB with W in the previous discussion. This is true 
for any climate diagnostic. More specifically, in the case of Tg, for W— j-SB 
and SB— T-W transitions. 



T,{S*„[C02]o)>T,{S},[C02 



J/J 



respectively. Note that for S* > 
transition to SB state can occur. 



and T,{S*, [C02]o) < T,{S}, [CO2]/) 

(12) 

1440 Wm ^, even if CO2 is ppm, no 



Figures 3(a) and 3(b) show the reconstruction for the materi al entropy 
.production, S^at p CS space, computed directly as described in Fraedrich 



and Lunkeit 



(2008 ). As with temperature, Smat increases monotonically with 



increasing S* and [CO2] on both manifolds. In the SB state, the entropy is 
mostly generated by dissipation of kinetic energy and by irreversible sensible 
heat transport, because the planet is almost entirely dry. For the W manifold 
the main contribution to entropy production comes from latent heat due to 
large scale and convective precipitation . In the bistable region, the range of 
Smat is (10, 19) lO-^W m-2 and (34, 62) lO'^W m'^K"! for the SB and 



W respectively, therefore a factor of 3 larger in the W regime with respect 
to the SB regime. This confirms that Smat may be a better indicator than 
temperature fo r discrim inating the SB and W states as already discussed in 
Lucarini et al. (2010b). Again, there is a range of values of Smat ^ from 19 
to 34 10~^W m~^ - which is not allowed by the system. 

In the bistable region, the SB and W states are quantitatively very differ- 
ent with respect to their physical properties. The two disjoint attractors can 
be thus thought to represent two different worlds, with completely different 
dynamical and thermodynamical properties. Therefore we treat them sepa- 
rately and then describe how the system makes a transition between them. 



15 



467 
468 
469 
470 
471 
472 
473 
474 
475 
476 



478 
479 
430 
481 
482 
483 
484 
485 
486 
487 
488 
489 
490 
491 
492 
493 
494 
495 
496 
497 
498 
499 
600 
501 
502 
503 



For this reason in the following two sub-sections dynamical and thermody- 
namical properties of the manifolds will be analysed individually in terms of 
the vertical and horizontal surface temperature differences, Carnot efficiency, 
meridional heat transport and dissipation of kinetic energy. Furthermore we 
shall relate these properties to the average mean global temperature and the 
material entropy production. As is illustrated as solid and dashed purple 
lines, in the CS space figures, each manifold will be divided up in to two sub 
regions: W, W/Bistable, on the upper manifold and SB, SB/Bistable on the 
lower. Then in a third section, we will analyse the transitions between the 
two manifolds occurring in bifurcations regions. 

4-2. The warm state 



The meridional energy transport profile is worked out, as explained in Lu- 



carini et al. (2011), by integrating over latitude the longitudinally averaged 
TOA radiation budget. We then define as a scalar indicator of the transport 
half of the sum of the peak values of the poleward heat transport in the two 
hemispheres, MET. In a moist atmosphere, the average global temperature 
and the meridional surface temperature difference, defined in our case as the 
difference between the mean surface temperature of the tropical (308, SON) 



and the polar (90S, 60S) and (60N,90N) regions (see Figure 5(a)), are the 
main contributing factors for controlling the meridional heat transport (see 



Figure 4(a)). This is due to the fact that temperature controls the latent 



heat released in the atmosphere because of the Clausius-Clapeyron effect 



iHeld and Soden, 2006) and the meridional temperature gradient controls 
baroclinicity of the atmosphere (Stone, 1978). Additionally, another mod- 



ulating factor is the vertical stratification of the atmosphere, as conditions 
of low stratification in the midlatitudes support stronger baroclinic activity 



for a given meridional temperature gradient (Holton, 2004). In the bistable 



region of the warm sector meridional heat transport has a fiat response to 
increasing S* and [CO2] and therefore T^. With increased T^, water vapour 
concentration of the atmosphere increases, thus leading to the strengthening 
of the poleward latent heat fluxes. In addition, the increased Tg causes sea 
and continental ice as well as seasonal snow cover to retreat towards the 
poles, thus lowering the surface albedo gradient. This contributes negatively 
to the changes in the meridional transport, because it leads to a decrease in 
the baroclinicity. 

In the W regime, the boundary between the bistable and the monostable 
regime approximately marks the point at which the Earth surface loses its 
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permanent sea-ice cover, thus supporting the idea that the presence of bista- 
bihty is intimately hnked with the powerful ice-albedo feedback. For Tg larger 
than approximately 300 K, the meridional temperature gradient decreases at 
a far slower rate with increasing Tg. This essentially means that in this re- 
gion the meridional heat transport is controlled only by the availability of 
water vapour in the atmosphere. Therefore at temperatures above 300K, the 
meridional heat transport becomes decoupled from the surface meridional 
temperature difference. This analysis consequently shows how important the 
involvement of the hydrological cycle is in the magnitude of the meridional 
heat transport and moreover, it indicates that when going from warm to very 
warm climates the hydrological cycle becomes the dominant climatic feature, 
leading to strong positive dependence of the meridional heat transport on the 



i^urface temperature. Our results agree with the findings of [Caballero and 



Langen (2005), who found in aqua planet simulations, that there is little 



518 scope for reducing the meridional temperature gradient further once the sea- 

519 ice and snow have melted. For a constant meridional temperature gradient, 

520 an increase in the meridional latent heat fluxes with increasing global mean 

521 temperature was shown. 

522 The conclusions drawn above find further support when looking at the 

523 mid-latitudes vertical temperature difference, ATy, defined as the mean tem- 

524 perature difference between the surface and the 500 hPa level (see Figure 



5(b) ) and which provides a rough measure of the stratification of the atmo- 



sphere of the atmosphere in the baroclinically active regime. The vertical 
temperature difference is largest along a band of the CS space centered half 
way through the bistable region in the direction of S*. Our current climate 
conditions would appear to be positioned at the centre of the band peak 
where such temperature difference peaks, implying conditions of reduced 
vertical stratification. For climates colder than present, increasing surface 
temperature causes the melting of sea-ice and of seasonal snow cover, so that 
the ensuing decrease in surface albedo (leading to increased surface absorp- 
tion) accounts for the increasing vertical temperature difference. Instead, in 
climates colder than present, the decrease in equatorial vertical temperature 
difference with Tg can be understood in terms of increased moist convection 
from warmer surface temperatures, resulting in an increase of moisture fluxes 
to the upper atmosphere, which then condense and release latent heat. In 
the bistable region, for the reasons discussed above, we expect to find a pro- 
nounced weakening of the dynamics of the climate system with increasing 
surface temperature. We test this hypothesis by computing the strength of 
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542 the Lorenz energy cycle (Fig. 4(b) ), which is equal to the average rate of dis- 

543 sipation of kinetic energy, and the Carnot efficiency of the system (Fig. |6(a)[ ), 

544 which measures, instead, how far from equilibrium the system is. The dis- 

545 sipation, similarly to the meridional temperature gradient, decreases mono- 

546 tonically with the surface temperature, and reaches its largest value just at 

547 the cold boundary of the W manifold. The efficiency is maximized along a 

548 band of climate states characterised by slightly higher temperatures. This 

549 region in the CS space lies about half way between the peak in the merid- 

550 ional and vertical temperature gradients, since a large value of efficiency, by 

551 its very definition, requires a compromises between maximizing these two 

552 gradients. The efficiency decreases monotonically with increasing CO2 or S*: 

553 warmer climates are characterised by smaller temperature differences, since 

554 the transport of water vapour acts as a very efficient means for homogenising 

555 the temperature across the system. The system has lower ability to produce 

556 mechanical work and is characterised by very strong irreversible processes, 

557 as described by the very large values of entropy production realised in these 

558 conditions. Consequently the value of the parameter of irreversibility a in- 

559 creases as conditions becomes warmer and warmer (not shown); see also 



560 



Lucarini et al. (2010b). 



561 4-3- The snowball state 

562 The SB state is intrinsically simpler than the W state because the hy- 

563 drological cycle has a negligible influence. This is due to the fact that at- 

564 mospheric temperatures are so low that we have in all cases an almost dry 

565 atmosphere. As a result we have, e.g. that, as opposed to the W state, 

566 the parameter a is small (~ 1) and weakly dependent on the value of the 



567 parameter S* and [CO2] (not shown); see also Lucarini et al. (2010b). 

568 Moreover, in the SB state the meridional gradients (and not only the 

569 globally averaged values) of albedo are very low and depend weakly on S* 

570 and the CO2, because under all conditions the sea surface is frozen almost 

571 everywhere and the continents are covered by ice and snow. The meridional 



572 heat transport (Figure 7(a) ) is much smaller than for the corresponding W 



573 states, because in the SB state a large fraction of the radiation is reflected 

574 back to space and the albedo gradients are small, so that the energy im- 

575 balance between low and high latitudes is small. Throughout the SB state, 

576 increases in S* and [CO2] , which lead to an increase in the surface temper- 

577 ature and in the rate of entropy production are accompanied by increases in 
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579 ure 
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581 



601 



608 



the meridional heat transport and in the dissipation of kinetic energy (Fig- 
7(b)[ ). This tells us that with increased meridional heat transport the 



intensity of the circulation also increases. However we find that the merid- 
ional temperature difference (not shown) has a very weak dependence on 

582 changing [CO2] and S*, as it varies by only 4 K across the explored param- 

583 eter range, so that the changes in the baroclinicity of the system cannot be 

584 due to changes in the meridional heat transport and in the intensity of the 

585 Lorenz energy cycle. The lack of large variations in the meridional surface 

586 temperature profile are essentially due to the fact that the net input of short- 

587 wave radiation is kept almost fixed by the constant surface albedo and by the 

588 virtue of almost no cloud cover. The system therefore can be seen as rather 

589 rigid, as changes in the absorbed radiation are almost exactly compensated 

590 by changes in the meridional heat transport. The re-equilibration mechanism 

591 must then contain elements which are not present in the classic baroclinic 

592 adjustment (Stone, [l978). We find that the midlatitudes vertical tempera- 



593 ture gradient (Figure 8(b) ) increases substantially with increased value of S* 

594 and [CO2], as an effect of the increased absorption of radiation near the sur- 

595 face. The reduced vertical stratification leads to more pronounced baroclinic 

596 activity even for a fixed meridional temperature gradient, thus leading to an 

597 increase in the meridional heat transport and in the intensity of the Lorenz 

598 energy cycle. The argument is in this case quite straightforward because we 

599 are considering an almost dry atmosphere. The efficiency has a dependence 

600 on S* and [CO2] which, as in the W case, is intermediate between those of 
the meridional and vertical temperature gradient fields. In the case of the 



602 SB state, the main signature is given by the vertical gradient (Figure 8(b)[ ) 



603 This further reinforces the idea that the investigation of meridional temper- 

604 ature gradients is not enough to grasp the mechanisms through which the 

605 system generates available potential energy and material entropy production 



606 ( jLucarini et al. 2011). 



4.4- Transition and comparison between manifolds 

Up to now, the W and SB states have been characterised as two entirely 

609 distinct climate regimes, and have underlined that the basic mechanisms 

610 of re-equilibration are rather different. In this subsection we would like to 

611 present some ideas aimed at making sense of the transitions between the 

612 two states occurring when we get close to the boundary of the upper or of 

613 the lower manifold. As has been seen, the W— J-SB transition is associated 

614 to a large decrease of surface temperature, rate of material entropy produc- 
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620 
621 



615 tion, and meridional heat transport. This is intimately related to the fact 

616 that whereas in the W state the hydrological cycle is a major player of the 

617 climate dynamics, in the SB state the hydrological cycle is almost absent. 

618 Nonetheless, this does not say much about the processes leading to the tran- 

619 sition between the two states, or better, describing how one of the attractors 
disappears. 

We have also discovered that the usual dynamical indicators of the at- 

622 mospheric state, i.e. the meridional temperature gradient and the vertical 

623 stratification do not necessarily indicate whether or not we are close to an 

624 irreversible transition of the system, e.g. by signaling something equivalent 

625 to a loss of elasticity of the system. In this regard, it is much more informa- 

626 tive to observe how the efficiency behaves near the transitions. We find that, 

627 as a general rule, each transition is associated to a notable decrease (more 

628 than 30%) of the efficiency of the system, and the closer the system gets to 

629 the transition in the CS space, the larger is the value of the efficiency. This 

630 can be interpreted as follows. If the system approaches a bifurcation point, 

631 its positive feedbacks become relatively stronger and the negative feedbacks, 

632 which act as re-equilibrating mechanisms, become less efficient. As a result, 

633 the differential heating driving the climate is damped less effectively, and 

634 the system is further from equilibrium, since larger temperature differences 

635 are present. Therefore, the system produces more work, thus featuring an 

636 enhanced Lorenz energy cycle and a stronger circulation. At the bifurcation 

637 point, the positive feedbacks prevail and the circulation, even if rather strong, 

638 is not able to cope with the destabilising processes, and the transition to the 

639 other manifold is realised. The new state is, by definition, more stable, and 

640 thus closer to equilibrium. The decreased value of the efficiency is exactly the 



642 



641 ipiarker of this property. This confirms what has been proposed by Lucarini 



et al.| ( |2010b ) but in much greater generality and can be conjectured to be 



643 a rather general properties of non-equilibrium systems featuring structural 

644 instabilities. 



645 5. Parameterisations 

646 In the previous sections, we have presented a systematic investigation 

647 of the thermodynamical properties of the planet in the SB and W states, 

648 obtained by exploring the CS parametric space. As mentioned above, in 

649 many cases the main quantity controlling the thermodynamic properties is 

650 the average surface temperature Tg only. This statement can be substantiated 
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651 by looking at Figs. [2]j8] and observing that in most cases the isohnes of 

652 the depicted thermodynamic quantities are approximately parallel, both in 

653 the W and SB state, to the isolines of Tg. This suggests the possibility 

654 of establishing approximate empirical laws of the form TssiS*, [CO2]) ~ 

655 TsBiTsiS*, [CO2])) and TwiS*, [CO2]) ^ TwiTsiS*, [CO2])), where T is a 

656 thermodynamical property such as, e.g. entropy production, and the lower 

657 index refers to whether we are in the SB or W state. For a given F, the 

658 empirical laws will be in general different in the SB vs in the W state. This 

659 result is quite interesting in a classical perspective of climate dynamics, where 

660 it is customary to parameterise large scale climate properties as a function of 

661 the surface temperature, especially when constructing simple yet meaningful 



678 



662 models (Saltzman, 2002) 



663 We want to complement this approach by testing whether it is possible 

664 to construct empirical laws expressing the thermodynamical properties of 

665 the system as a function of its emission temperature Te = {LWtoa/o')^^'^ 

666 - where LWtoa is the area- and time-averaged outgoing long wave radia- 

667 tion at the top of the atmosphere and a the Stefan-Boltzmann constant - , 

668 which is a quantity of more readily astrophysical interest and more funda- 

669 mental physical significance. Therefore, we would like to be able to observe 

670 well-defined functions of the form Tsb{S% [CO2]) ^ Tsb{Te{S*, [CO2])) and 

671 Tw{S*, [CO2]) ~ Tw(Te{S*, [CO2])), which, given the observation made 

672 above, requires also the existence of an approximate relationship Te ~ Te(Ts). 

673 Indeed, Fig. [9]shows that it is possible to observe an empirical, monotonic 

674 relationship between and Te in a vast range of climate states, including 

675 both the W and SB states. Obviously, the correspondence between the two 

676 variables is not perfect, so that for each value of Tg it is possible to have 

677 climates featuring a range of values of Te, and vice versa. The spread asso- 
ciated with the ensemble points is of order 10 K and 5 K for Te, and 20 K 



679 and 10 K for Ts in the W and SB states respectively. A range of values of 

680 Te (230 — 240 K) and Tg (240 — 260 K) are not permitted by the system as 

681 a result of the ice-albedo instability. We can now proceed to see to which 

682 extent the quantities and Te can be used as a good predictor of the ther- 

683 modynamic properties of the planet and compare their skill with Tg. The 

684 following description is, on purpose, purely qualitative, because we mostly 

685 want to show to what extent Te and Tg can be used as predictor, rather than 

686 discussing in detail what are the appropriate transfer functions. 



687 



Figure 10(a) and 10(b) show that a clear monotonic relationship exists 
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between Tg, Tg and the material entropy production Smat, with higher val- 
ues of Te and Tg associated with larger values of Smat- It is important to 
note that while the values of Smat and Tg are related to a very high degree 
of approximation on a one-to-one basis, the functional dependence between 
Te and Smat is less well defined (the spread of points is much wider than in 
the Tg case, and especially so in the very warm regimes). Note also that the 
dependences, while monotonic in both cases, are very different between the 
W and SB states. The common line is that the more is the energy absorbed 
by the planet (larger values of Te), the more intense in the production of 
material entropy. Note that this is indeed not the case when considering the 
intensity of the Lorenz energy cycle (Fig. 11(a) and ll(b)[ ): higher absorp- 
tion of energy corresponds to a weaker energy cycle in the W state and to 
a stronger cycle in the SB state. Therefore, the increase in the entropy pro- 
duction with Te is due, in the SB case, to the increased dissipation of kinetic 
energy, while in the W state it is due to the greatly enhanced hydrological 
cycle. In agreement with what has been discussed in the previous section, one 
can explain the contrasting responses of the W and SB states to increases in 
the absorption of energy by analyzing the Carnot efficiency of the planetary 



fluid (Fig. 12). It is interesting to note that, at the transitions, the values of 
the efficiency are practically equal: rjw^sB ~ Vsb,w ~ 0.03. Moreover, we 
observe that t] saturates in the very cold regime of SB states and very warm 
regime of W (see Fig. 12) at rjgat^w ~ Vsat,SB ~ 0.012. Also for the Lorenz 



energy cycle and the efficiency, T, does a much better job as predictor of the 
thermodynamic quantities, since the spread of points is much higher when 
considering Te- 

Last, we want to present the link between Te, Tg and the meridional 
energy transport MET (Fig. 13). We observe that in the SB regime there 
is a weakly positive relationship between Te (Tg) and the transport, for the 
reasons described in the previous section; instead, for a vast range of values 
of Te (Tg) in the W regime, the transport is almost insensitive to Te (Tg), 



discu ssed thoroughly in Caballero and Langen (2005) and Donohoe and 



Battisti (2011 2012). The rigidity of the climate system has been attributed 



by Donohoe and Battisti (2011, 2012) to the fact that the meridional heat 



transport is mainly determined by planetary albedo and thus atmospheric 
composition rather that surface albedo and sea-ice coverage. However, such 
a mechanism ceases to exist when the ice albedo feedback becomes ineffective 
because of the disappearance of sea ice above a threshold value of Te ~ 255 
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742 
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725 K {Tg ^ 300 K). Above this value, we observe a steep monotonic increase of 

726 the transport with temperature, because changes in the latent heat transport 

727 are mainly responsible for this behavior. This agrees with the idea that the 

728 dynamics and sensitivity of a warm planet is, in some sense, dominated by 

729 the hydrological cycle. In this latter case, the skill of the two temperature 

730 quantities in parameterising the MHT is comparable. 

731 Figs [T0HT3I clearly define the transfer functions allowing to compute the 

732 thermodynamic properties of the planet from the value of Te or Tg only: 

733 this is definitely intriguing because it paves the way for computing non- 
734 equilibrium properties from a quantity describing the zero-order radiation 

735 balance of the system at the top of the atmosphere and of the radiative- 

736 convective equilibrium at surface (in the case of T<j). Let us finally note that, 

737 in general, the functional dependences of the thermodynamic quantities on 

738 Tg is more precisely defined than the one on Te since these show a larger 

739 spread of the ensemble points. This is particularly evident in quantities such 
as Smat (Fig. 10) and ?7 (Fig. 1^, W (Fig. n)) but less in others, e.g. MET 



741 (Fig. 13). Figures (|9][12|) leave various open questions which definitely need 
to be accurately addressed in future studies: i) why is it possible to use 
to a good degree of approximation the quantity Te to parameterise the non- 
744 equilibrium thermodynamical properties of a system? ii) why does the surface 

745 temperature does an even better job (actually, astounding) in controlling such 

746 complex properties? iii) is the result on the decrease of the efficiency at the 
bifurcation points in both W— )-SB and SB— i-W transitions of more general 



relevance? 



749 6. Summary and conclusions 

750 Motivated by the recent discoveries of extrasolar planets, in this contribu- 

751 tion we have shown that the habitability conditions, for an Earth-like planet, 

752 is a necessary but not solely sufficient condition since two stable sets of states 

753 are possible, with one of the two characterized by the global glaciation of wa- 

754 ter. In particular we have shown how the stability properties of the climate 

755 system depend on the modulation of the two main parameters describing 

756 the radiative forcing, i.e. the energy input of the parent star and the long 

757 wave opacity of the atmosphere. In our analysis we propose that the point of 

758 view of non-equilibrium thermodynamics is especially useful for understand- 

759 ing the global properties of the climate system and for interpreting its global 

760 instabilities. 
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We have discovered that in a rather wide parametric space, which in- 
cludes the present chmate conditions, the chmate is multistable, i.e. there 
are two coexisting attractors, one characterised by warm conditions, where 
the presence of sea ice and seasonal snow cover is limited (W state), and 
one characterised by a completely frozen sea surface, the so-called snowball 
(SB) state. Our results, obtained using the climate model PlaSim, confirm 
and extend what has been obtained in various studies performed with mod- 
ipls of various degrees of complexity. We point the reader to [Pierrehumbert 



et al. (2011) and Lucarini et al. (2010b) for an extensive discussion. In this 



regard, the main improvement of this work is that a two dimensional para- 
metric space is explored (whereas usually variations in the solar constant or 
in the opacity of the atmosphere are considered separately), which allows the 
gathering of more complete information on the possible states of the climate 
system and actual mechanisms relevant in a paleoclimatological context, as 
explained in Pierrehumbert et al. (2011). For all considered values of [CO2], 
which range from 90 to about 3000 ppm, the width of the bistable region is 
about 200 Wm~^ in terms of the value of the solar constant, and its posi- 
tion depends linearly on the logarithm of the [CO2], being centered around 
smaller values of the solar constant for increasing opacity of the atmosphere, 
shifting by about 15 Wm~^ per doubling of CO2 concentration. The W state 
is characterized by surface temperature being 40 K-60 K higher than in the 
SB state, and also, in terms of the material entropy production, is larger by a 
factor of 4 (order of 40 - 60 lO'^Wm'^ vs. 10 - 15 lO'^Wm'^). The bound- 
aries of the bistable region are approximately isolines of the globally averaged 
surface temperature, and in particular, the warm boundary, beyond which 
the SB state cannot be realized, is characterized by vanishing permanent 
sea ice cover in the W regime. This reinforces the idea that the ice-albedo 
feedback is the dominant mechanism for the multistability properties. 

The thermodynamical and dynamical properties of the two states are 
largely different, as if we were discussing two entirely different planets. In 
the W state the climate is dominated by the hydrological cycle and latent 
heat fluxes are prominent in terms of redistributing the energy in the system 
and as contributors to the material entropy production. The SB state is 
eminently a dry climate, where heat transport is realized through sensible 
heat fluxes and entropy is mostly generated through the dissipation of the 
kinetic energy. The dryness of the SB atmospheres also explains why the 
climate sensitivity is much smaller. 

In the bistability region of the W states, the meridional heat transport 
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is rather constant, as the contrasting effect of the enhancement of latent 
heat fluxes driven by increasing surface temperature and the reduction in 
the baroclinicity due to the decrease in the meridional temperature gradient 
compensate almost exactly. In the W state beyond the bistability region, 
the meridional heat flux increases quite dramatically with the surface tem- 
perature (Fig. 5(a) and 13) since the compensating albedo mechanisms are 
shut off as sea ice is completely removed from the surface. In the SB states, 
increased incoming radiation or increased [CO2] leads to increases in the 
meridional heat transport. In this case, the water vapor plays little role, 
and, somewhat surprisingly, the meridional temperature difference has also 
a rather flat response to the parametric modulation. In this case, the dom- 
inant mechanism determining the properties of the meridional heat transfer 
is the change in the vertical stratification, which becomes weaker for warmer 
climate conditions. This implies that the atmospheric circulation strength- 
ens for increasing values of the solar constant and of [CO2]. In fact, the 
Lorenz energy cycle becomes stronger for warmer climate conditions, and 
the Carnot-like efficiency of the climate system has an analogous behavior. 
The opposite holds for the W state, where the intensity of the Lorenz energy 
cycle and the efficiency decrease for warmer conditions, the reason being that 
the water vapor becomes more and more efficient in homogenizing the system 
and destroying its ability to generate available potential energy. 

A general property which has been found is that, in both regimes, the 
efficiency increases when we get closer to the bifurcation point and at the 
bifurcation point the transition to the newly realized stationary state is ac- 
companied by a large decrease in the efficiency. This can be framed in a rather 
general thermodynamical context: the efficiency gives a measure of how far 
from equilibrium the system is. The negative feedbacks tend to counteract 
the differential heating due to the sun's insolation pattern, thus leading the 
system closer to equilibrium. At the bifurcation point, the negative feedbacks 
are overcame by the positive feedbacks, so that the system makes a global 
transition to a new state, where, in turn, the negative feedbacks are more 
efficient in stabilizing the system. 

The results discussed in this paper support the adoption of new diag- 
nostic tools based on non-equilibrium thermodynamics for analysing the 
fundamental properties of planetary atmospheres. The next step in this 
direction is a more quantitative understanding of the global relationships 
between surface temperature, material entropy production, meridional heat 
fluxes and Carnot-like efficiency for the SB and W states and to propose pos- 
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837 sible parametrisations for the dynamical and thermodynamical properties of 

838 the SB and W states. Another hne of research will explore the dependence 

839 of these quantities on other fundamental parameters, e.g. the rotation rate 

840 and the surface drag, relevant for atmospheres of terrestrial planets. Re- 

841 cent theoretical results suggest that it is possible to derive up to a good 

842 degree of approximation the full thermodynamical properties of planetary 

843 atmospheres from the coarse observations of radiative fluxes at TOA and 



844 surface (Lucarini et al. , 2010b). In this work, we have pursued this line of 

845 research by proposing well-defined empirical functions allowing for express- 

846 ing the global non-equilibrium thermodynamical properties of the system in 

847 terms of the globally averaged surface temperature and planetary emission 

848 temperature only. The skills of such temperature quantities as predictors 

849 is fairly good. The reasons why these parameterisations work so well, and 

850 why this is especially the case for globally averaged surface temperature def- 

851 initely require further investigations. These findings fully confirm the results 

852 discussed in Lucarini et al. ( 2010b|[a ) and generalize them to a much larger 



853 variety of climates, explored by changing both the energy input of parent 

854 star and the opacity of the atmosphere. These results reinforce the relevance 

855 of the globally averaged surface temperature as a meaningful thermodynam- 

856 ical quantity of the climate system, beyond its obvious practical importance. 

857 Moreover, they pave the way for the possibility of practically deducing fun- 

858 damental properties of planets in the habitable zone from a relatively simple 

859 observable. 
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Figures' captions 
• Figure [T] 

(a) Surface temperature and (b) material entropy production for steady 
states obtained for different values of the solar constant S*. The present 
climate is marked with a black circle, the W states in red and the SB 
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states in blue. Adapted from Lucarini et al. (2010b) 



Figure |2] 

Contour plot of surface temperature (K) as a function of S* and the 
[CO2]. The lower SB (a) and upper W (b) manifolds are shown. The 
transition SB— )-W and W— )-SB are shown by the upper and lower purple 
lines respectively. The blue dots indicate the values of (S'*,[C02]) for 
which simulations have been performed. 

Figure |3] 

Contour plot of Smat (10"^Wm"^K"^) as a function of S* and the [CO2] 
for the lower SB (a) and upper W (b) manifolds. 



1004 • Figure |4] 

1005 Contour plot of: (a) meridional energy transport MET (PW, 1 PW= 
10^^ W) and (b) average rate of dissipation of kinetic energy W (W 
m~^) as a function of S* and the [CO2 ] for the W states. 



• Figure |5] 

Contour plot of: (a) meridional temperature gradient (K) and (b) mid- 
latitude vertical temperature difference (K) as a function of S* and the 

1011 [CO2 ] for the W states. 

1012 • Figure |6] 

1013 Contour plot of the efficiency as a function of S* and the CO2 concen- 

1014 tration for the W states. 

1015 • Figure [7] 

1016 Contour plot of: (a) meridional energy transport MET (PW) and (b) 

1017 average rate of dissipation of kinetic energy (W m~^) as a function of 

1018 S* and the [CO2 ] for the SB states 

1019 • Figure |8] 

1020 Contour plot of (a) the Carnot efficiency rj and (b) midlatitude vertical 
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temperature difference (K) as a function of S* and the [CO2 ] for the 
SB. 

Figure |9] 

Surface temperature (K) versus emission temperature Te (K) for the 
W and SB states (each circle represents a simulation). We have marked 
the transitions W — )■ SB (red dot - arrow - blue dot) and SB — )■ W 



(blue dot - arrow - red dot). Same convention is used in Figs. 10 13 



Figure 10 



Material entropy production Smat (10 m ^ K ^) vs. (a) emission 
temperature Te (K) and (b) surface temperature Tg (K). 



Figure 11 



Lorenz Energy cycle strength W (W m ^) vs. (a) emission temperature 
Te (K) and (b) surface temperature (K). 



Figure [12] 

Carnot efficienct 77 vs. (a) emission temperature Te (K) and (b) surface 
temperature Ts{K). 



Figure 13 



Meridional heat transport index MET (in PW) vs. (a) emission tem- 
perature Te (K) and (b) surface temperature Tg (K). 
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